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q ' Abstract. YBa2Cu30y 24° (30°) bicrystal grain boundary junctions (GBJs), 

O . shunted with 60 nm (20 nm) thick Au, were fabricated by focused ion beam milling 

with widths 80 nm < w < 7.8 /im. At 4.2 K we find critical current densities j c in 
the 10 5 A/cm 2 range (without a clear dependence on w) and an increase in resistance 
, times junction area p with an approximate scaling p cx w 1 / 2 . For the narrowest GBJs 

jy-^ ' j c p w 100 pV, which is promising for the realization of sensitive nanoSQUIDs for the 

\& • detection of small spin systems. We demonstrate that our fabrication process allows 

the realization of sensitive nanoscale dc SQUIDs; for a SQUID with w w 100 nm wide 
0^ ' GBJs we hnd an rms magnetic flux noise spectral density of S^ 2 « 4/i$o/Hz 1 ^ 2 

Oin the white noise limit. We also derive an expression for the spin sensitivity S*// 2 , 
1 /2 

which depends on , on the location and orientation of the magnetic moment of a 
magnetic particle to be detected by the SQUID, and on the SQUID geometry. For the 

,1/2 _ 

S3 



not optimized SQUIDs presented here, we estimate = 390 pb/ V Hz, which could 
be further improved by at least an order of magnitude. 



PACS numbers: 85.25. CP, 85.25.Dq, 74.78.Na, 74.72.-h 
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1. Introduction 

There is a growing interest in developing sensitive miniaturized superconducting 
quantum interference devices (SQUIDs) for the investigation of small spin systems and 
scanning SQUID microscopy with sub- /mi spatial resolution [HE! [3]. A main motivation 
are measurements on single nanomagnetic particles, and the ultimate goal is the direct 
detection of switching of a single electronic spin with various potential applications 
in spintronics, quantum computing and on biomolecules. These applications require 
sub-/im Josephson junctions and SQUID loops, both for optimum inductive coupling to 
nanosized objects and for improving the SQUID sensitivity for the operation at switching 
fields of the magnetic particles up to the Tesla range at temperature T « 4 K and below 
(see e.g. Refs. [21 HJ El Ej and references therein). 

Most frequently used techniques for creating superconducting sub-/im thin film 
structures are based on electron beam lithography, or on focused ion beam (FIB) 
patterning. While the patterning of sub-/im SQUID loops poses no particular problems, 
the properties of sub-/zm Josephson junctions embedded in the SQUID loop have to be 
carefully optimized in order to realize low-noise SQUIDs. Here, an important figure of 

1/2 

merit is the rms spectral density of flux noise S$ , which for optimized SQUIDs based on 
the standard Nb/AlO^/Nb technology for /xm-sized junctions is in the range of /x$o/ VTIz 
($o is the magnetic flux quantum). However, critical current densities j c > lkA/cm 2 
are hard to achieve reliably with this technology, yielding too small critical currents 
I c = j c Aj for junction areas Aj in the (100 nm) 2 range. Therefore, recent research 
focused on Nb or Al thin film constriction type junctions (with widths u>£50nm), 
sometimes shunted with a thin film normal metal layer (e. g. Au, or W) to ensure non- 
hysteretic current- voltage-characteristics (IVCs) [3J. This approach typically produced 

1 /2 / 

SQUIDs with S% *. 1 m$ /v Hz. Only very recently, highly sensitive Nb SQUIDs with 
S$ ~ 0.2/i$ /\/Hz have been demonstrated [7]. These encouraging results have been 
obtained with FIB patterned constriction junctions at T = 6.8 K. 

For operation in high magnetic fields B a small junction size (_L B) is required due 
to the suppression of I c with increasing magnetic flux in the junction above several $o> 
which demands for a particularly large j c . Furthermore, the maximum field of operation 
is limited by the upper critical field B C 2, e. g. to £ 1 T for typical Nb thin film based 
SQUIDs. Here, high-transition-temperature (high-T c ) SQUIDs offer three advantages: 
(i) very high B c i in the tens of Tesla regime or even more, (ii) high j c > 10 5 A/cm 2 
for grain boundary junctions (GBJs) operating at T = 4.2 K and well below and 
(iii) a GBJ geometry with the junction barrier perpendicular to the thin film SQUID 
loop; this allows to apply very large in-plane fields (for switching the magnetization of 
nanoparticles) which do not couple to the SQUID and which do not reduce I c . Here, 
the challenge is to produce sub-/im GBJs with high quality, in particular with high j c 
and p = RAj, where R is the junction resistance. At this point, we should mention, 
that in principle, also constriction type junctions based on HTS, or MgB2 thin films 
may fulfill the above mentioned requirements, and in fact such junctions and SQUIDs 
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based on them have been patterned by FIB [HI El EU Ell EE2]- However, these junctions 
typically often show flux-flow type or hysteretic IVCs at low T, such that their operation 
temperature is often limited to a narrow T range well above 4K, and their performance 
with respect to flux noise so far did never reach the performance of high-T c low-noise 
GBJ SQUIDs [B]. 

Already in the 1990s thin film high-T c YBa 2 Cu 3 7 (YBCO) sub-/im GBJs and 
SQUIDs have been fabricated using e-beam lithography [TH [15]. However, oxygen 
loss during processing, in particular for very thin films, may require post-deposition 
annealing to improve junction characteristics [16]. More recently, sub-/xm YBCO GBJs 
have also been fabricated by FIB [17] . and both technologies enabled fundamental 
studies on transport and noise in high-T c sub-^m GBJs [161 HH1 HH1 ED] . Still, a significant 
degradation of j c for w S 500 nm was found [HI HSl [17] , and the use of deep sub-/xm GBJs 
for the realization of nanoSQUIDs has not been explored yet. The motivation for the 
realization of sub-/im GBJs with widths well below 500 nm is based on the following 
considerations: First of all, operation in high magnetic fields in the Tesla range as 
mentioned above, requires very accurate alignment of the applied magnetic field in the 
thin film plane in order to reduce coupling of the applied out-of-field component to the 
GBJ. This demands for as small as possible GBJ widths. Furthermore, as will be shown 

1/2 

below, the spin sensitivity scales linearly with the rms flux noise S£* of the SQUID. 
Optimization of S$ requires an as small as possible SQUID loop inductance L [21J, 
i.e. minimization of the dimensions of the SQUID loop, which for topological reasons has 
to be intersected by the grain boundary. Here, SQUID loop sizes of the order of 100 nm 
seem to be feasible, according to our experience on FIB patterning of our devices as 
described below. In order to ensure optimum SQUID performance, one should achieve 
at least a few SQUID modulations within the Fraunhofer-like I C {B) modulation of the 
single GBJs. This in turn requires also shrinking the GBJ widths down to the size of the 
SQUID loop. Hence, our goal is to demonstrate the feasibility of FIB patterning YBCO 
GBJs down to junction widths of the order of 100 nm. In order to accomplish this, we 
investigated the scaling behavior of YBCO GBJ properties with linewidths ranging over 
two orders of magnitude, from ~ 8 /im down to 80 nm, and we investigated the electric 
transport and noise properties of YBCO GBJ dc SQUIDs with the smallest linewidths 
achieved within this study. We note that we performed so far only investigations on low- 
field properties of the fabricated GBJs and SQUIDs, as we are at this stage interested 
in clarifying the intrinsic scaling properties of our devices with GBJ width, although 
the ultimate goal of this work is to operate such SQUIDs in high magnetic fields in the 
Tesla range at T = 4.2 K and well below. 

The remainder of this paper is organized as follows. Section 2 very briefly addresses 
sample fabrication and layout, including some information on the quality of our YBCO 
thin films. In Sect. 3 we first describe and discuss the results of electric transport 
properties of our shunted GBJs, with focus on their dependence on junction width, 
which was varied over two orders of magnitude (3.1). The second part (3.2) of this 
section describes the results obtained for our SQUIDs, with focus on electric transport 
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and noise properties of the SQUID with the smallest GBJ width used in this study 
(SQUID 2). Having characterized our SQUIDs, we discuss in Sect. 4 the important 
relation between the flux noise S^ 2 of the SQUIDs and the spin sensitivity S l J 2 , which 
is the important figure of merit for detection of small spin particles. Here, we provide a 
solution for calculating the spin sensitivity for any arbitrary geometry of the SQUID loop 
as a function of position and orientation of the magnetic moment of a small particle 
to be detected. We then apply this solution to the particular geometry of SQUID 2 
and finally discuss perspectives for further optimization of S l J 2 . Section 5 contains our 
conclusions. 

2. Sample Fabrication 

We fabricated devices on SrTiOs (STO) symmetric [001] tilt bicrystal substrates with 
misorientation angle 8 = 30° (chip 1) and 24° (chip 2). Figure [U(a)-(c) illustrates 
the fabrication steps. We deposited d Y = 50 nm thick c-axis oriented epitaxially grown 
YBCO by pulsed laser deposition (PLD), followed by in-situ evaporation of Au with 
thickness d^ u — 20 nm (chip 1) and 60 nm (chip 2), serving as a resistive shunt and 
protection layer during FIB milling. For details on PLD growth of our YBCO films 
on STO substrates, and their structural and electric transport properties see Ref. [22] . 
In brief, our 50 nm thick YBCO films typically yield 0.1° full width half maximum of 
the rocking curve at the (005) x-ray diffraction peak, have T c = 91 K with a transition 
width ~ 0.5 K and resistivity p ~ 50 pQcm at T = 100 K. On both chips, 8 fim wide 
bridges straddling the grain boundary were fabricated by photolithography and Ar ion 




GB GB GB 



Figure 1. (Color online) Upper row: Schematic illustration of the steps used for 
fabricating YBCO grain boundary junctions (GBJs). (a) in-situ deposition of a 
YBCO/Au bilayer on a bicrystal STO substrate; (b) patterning of 8 //m wide bridges 
straddling the grain boundary (GB) (by photolithography and Ar ion milling); (c) 
patterning of a narrow GBJ by FIB. The location of the GB is indicated by dashed 
lines. Bottom row (d)-(f) shows scanning electron microscopy images of three single 
junction devices (from chip 1) with different width (indicated in the graphs). 
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Figure 2. SEM images of the three SQUIDs (loop size 1.0 x 1.2 /im 2 ) fabricated on 
chip 2. Labels in black boxes give junction widths. 



milling and then patterned by FIB with Ga ions (50 pA, 30 kV) to make junctions and dc 
SQUIDs with junction widths 80 nm < w < 7.8 fim. FIB patterning was performed with 
a dual beam 1540 XB cross beam (Zeiss). This allowed us to apply an optimized FIB 
cut procedure (soft FIB procedure), with small ion current density and minimum ion 
exposure time of non-milled areas, i.e. only very brief snap shot imaging prior to milling. 
Even for imaging by the electron beam, we minimized the exposure time in order to avoid 
damaging of our FIB cut bridges. Figure [TJd)-(i) shows scanning electron microscopy 
(SEM) images of three GBJs fabricated on chip 1. In total, we investigated 22 single 
GBJs and 3 dc SQUIDs (on chip2; hole size 1.0 x 1.2 /im 2 ). SEM images of the three 
SQUIDs are shown in Fig. [2j 



3. Experiments 

We characterized our devices at T = 4.2 K in a magnetically shielded environment. For 
measurements of IVCs, I C {B) and V(B) we used a 4-point arrangement with a room 
temperature voltage amplifier. For SQUID noise measurements we preamplified the 
output signal with a Nb dc SQUID amplifier with O.lnV/vTIz resolution. 



3.1. Transport properties vs junction width 

All devices showed resistively-and-capacitively-shunted-junction (RCSJ)-type IVCs, 
which for some of the sub-pm junctions on chip 1 (thinner Au shunt) had a small 
hysteresis. Therefore, for chip 2, we increased d\ u by a factor of three, yielding non- 
hysteretic IVCs, except for the 530 nm wide junction, which has an exceptionally high 
j c p. From the IVCs we determined I c , and R, and calculated j c and p, using dy = 50 nm 
and w as obtained from SEM images. The results of these measurements are summarized 
in Fig. El plotted vs w which spans two orders of magnitude. Full symbols in Fig. [3] 
show data obtained after FIB patterning; open symbols show data (for chip 2) from the 
8 /jm wide bridges prior to FIB patterning. 

Figure [3]^a) shows j c {w), which is well above 10 5 A/cm 2 over the entire range of 
w. The j c values shown are typical for 9 = 24° and 30° YBCO GBJs at 4.2 K and 
w /im[T3]; however, such high j c values have not been previously observed for widths 
down to 80 nm. We do find a significant scattering of j c (w), however without a clear 
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Figure 3. (Color online) Transport data of YBCO GBJs and SQUIDs vs junction 
width w (solid symbols): (a) critical current density j c {w); (b) junction resistance 
times area p(w); (c) j c p(w). Dashed (solid) lines indicate average j c values [in (a)] and 
approximate scaling of p(w) [in (b)] and j c p(w) [in (c)] for chip 1 (2). Open squares are 
data for the same GBJs on chip 2, measured prior to FIB patterning (i.e. w = 8 pm). 
The numbers on the top axis in (a) label the device numbers on chip 2. 



width dependence. The average j c for chip 2 is 1.5 times the one for chip 1, as expected 
from the scaling j c {&) of GBJs [23]. The comparison of j c of the same bridges before 
and after FIB patterning shows that for most devices j c even slightly increased after 
FIB patterning. The position of the devices on chip 2 (along the GB of the substrate) 
is ordered according to their device number (1-15) [c.f. top axis of Fig. 13(a)] from the 
left to the right edge of the substrate. There is a clear trend of increasing j c by about a 
factor of two (for the 8 /im wide bridges) along the entire substrate. The origin of this 
gradient in j c has not been clarified; however we can rule out a corresponding variation in 
the YBCO film thickness across the substrate. A possible explanation for the observed 
gradient in j c of the 8 /im wide GBJs along chip 2, could be a gradient in the quality 
of the GB in the bicrystal substrate, which in turn can cause a gradient in the barrier 
thickness of the GBJs along the chip. 

Figure E(b) shows an approximate scaling p oc y/w of unclear origin. We note that 
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the lines shown in Fig. E^b) and Fig. [3^c) are no fits to our data. These lines are just 
drawn to illustrate the trend of decreasing resistivity p with decreasing junction with w. 
Before FIB patterning, the 8 pm wide GBJs on chip 2 all had p 0.17f2/im 2 , which falls 
onto the observed p{w) dependence, indicating that this scaling is not specific to FIB 
patterned GBJs. Furthermore, p ~ 0.17f2/im 2 is an order of magnitude below typical 
values for unshunted GBJs [23], which we attribute to the Au shunt, and which is also 
consistent with the larger p of GBJs on chip 1 with thinner Au. For chip 1, d^u = 20 nm 
is close to the 15 nm implantation depth of 30 keV Ga ions in Au [21] . Hence one might 
expect that FIB induces an increase in the Au resistivity via Ga implantation. This 
effect should be suppressed for chip 2 with 3 times thicker Au. In any case, it is hard 
to explain, why Ga implantation should increase p for wider junctions. Certainly, Ga 
implantation is not the only detrimental effect of FIB patterning. In particular, the Ga 
beam might destroy the crystalline order close to the patterned edges. However, our 
experimental observation of almost constant j c for GB J widths down to 80 nm rules out 
severe edge damage effects on a length scale of several tens of nm. This observation 
also rules out such effects as a possible explanation for the observed scaling behavior of 
p{w). 

Figure [3](c) shows j c p(w) ~ 0.1 . . . lmV, i.e. at least one order of magnitude below 
the values for unshunted YBCO GBJs of comparable j c . This is certainly due to the 
suppression of j c p by the Au shunt required to ensure non-hysteretic IVCs at 4.2 K. The 
decrease of j c p with decreasing w is due to the scaling of p{w) mentioned above. Still, 
even for the 80 nm wide GBJs we find reasonable values of j c p around 100 /iV, which 
are certainly quite suitable for the realization of sensitive SQUIDs. 

3.2. SQUID parameters, transport characteristics and noise performance 

The results of transport measurements on all three SQUIDs (on chip 2) are summarized 
in Tab. [D The SQUID inductance L was calculated with the numerical simulation 
software 3D-MLSI [25], which is based on a finite element method to solve the London 
equations for a given film thickness and London penetration depth {dy = 50 nm and 
\l = 140 nm, respectively, in our case). As dy <C A^, the kinetic inductance contributes 
significantly to L. For SQUID 1 and SQUID 2, the GBJ widths Wi (i = 1,2) are 
below Xl, which increases L over the one of SQUID 3 with wider junctions. From the 
calculated L and measured I c we obtain (3l = LJ c /$o ~ 0.3. . .0.7, i. e. not far from 
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Table 1. Parameters of YBCO GBJ dc SQUIDs. 
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Figure 4. (Color online) Transport and noise characteristics of SQUID 2: (a) IVC at 
B = 0; inset shows I C (B). (b) V(B) for / = -70.8 ... - 12.3 and 10.0 . . . 68.5 fiA (in 
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3.9 /iA steps), (c) Spectral density of mis flux noise S$ (/); inset: SEM image of the 
SQUID, (d) V(<$>) and S$ /2 ($) (averaged from / = 5 to 6kHz). 



the optimum value (3l ~ 1 [26]. The transfer function V$, i. e. the slope of the U($) 
curves at optimum bias current and applied flux $ = ±^$o is around 0.1mV/$o, and 
the effective area A e fj = <&/B « 8 /mi 2 for all three SQUIDs. 

Figure |H shows electric transport and noise data obtained for SQUID 2 [the device 
with smallest w; see inset in (c)]. Figure H^a) shows an IVC for an applied field B = 
corresponding to a maximum in I c . The small jump at I c to V ^ indicates that 
the junctions are at the transition to the underdamped regime. The inset in Fig. 11(a) 
shows I C {B) with 40% modulation. Figure IU(b) shows V(B) for various bias currents 
I. The small shift in the minima of V(B) upon reversing I is in accordance with the J c 
asymmetry of the two GBJs due to their different widths. 

Finally, graphs (c) and (d) in Fig. H] show the results of noise measurements on 

1 /2 

SQUID 2. Fig. H](c) shows the rms spectral density of flux noise S$ (/) oc f~ x for 
optimum flux bias $ = — 0.286 $o with x ~ 1.8 for frequencies /^lkHz, which we 
attribute to J c fluctuations in the GBJs p3]. For larger / we find a white flux noise 

1/2 / 1/2 

level ^ ^ ~ 4 /i$o/ v Hz, which to our knowledge is the lowest value of S$ obtained for 

TTV 1 /2 

a YBCO dc SQUID with sub-/im GBJs so far. Fig. HJd) shows the rms flux noise S$ 
(averaged from / = 5 to 6 kHz) and the SQUID voltage V vs applied flux $. We find 

1 /2 

a rather shallow minimum in S$ ($) for an applied flux where the slope of the V($) 
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curve [also shown in graph (d)] is close to its maximum. 
4. Spin sensitivity 

Coming back to the main motivation of this work, i.e. the development of nanoSQUIDs 
for the detection of small spin systems, we derive an expression for the spin sensitivity 
SjJ 2 , which we then use to calculate S l J 2 for the particular geometry and flux noise of 
SQUID 2 as a function of the position of a magnetic particle for a given orientation of 
its magnetic moment. is the spectral density of spin noise, which depends on the 
spectral density of flux noise S$ of the SQUID and on the coupling between a magnetic 
particle with magnetic moment ft — ji ■ and the SQUID via the relation = S^/cj) 2 ^. 
Here, (^(e^, r M ) = ^(/Z, i s the magnetic flux $ M per magnetic moment /i coupled 

into the SQUID loop by the magnetic particle, which is located at the position r M and 
which is oriented along e M . This means, in order to determine S l J 2 for a given S^ 2 , 
one needs to calculate the coupling function (f)^(e^, r^), which will also depend on the 
SQUID geometry. 

To determine (j)^, we assume that the magnetic moment /7 is moved from a distance 
far away to a position r = r M close to the SQUID loop. When the magnetic moment 
approaches r M , a circulating current r M ) is induced in the SQUID loop, which 
compensates the coupled flux due to the diamagnetic response of the SQUID loop. 
The magnetic field energy stored in the loop of inductance L is Wi oop = ^LI 2 . The work 
required to place the particle in the magnetic field B^if) produced by the circulating 
current 1^ is W M = — hjl- B^I^r^). We note that > 0, due to the diamagnetic 
response of the SQUID loop. Hence, the total work required to bring the magnetic 
particle to the position r M is 

w 1 = w loop + Wfj, = \ul - ■ rg . (i) 

On the other hand, instead of the SQUID, we may consider a fixed current system 
producing the same field B^il^^r^) as the SQUID, when the particle is in its final 
position ry In this case, the particle has a (positive) energy 

W2 = -/2-4(^V)- (2) 
From Wx = W 2 we obtain I 2 = -fl- B^/L. With ^ = LI^ and with B^/I^ = B/I = b 
one thus obtains 

^iBllA = M e„ f>) = -S„ • b(f,) , (3) 

where / is an arbitrary current circulating in the SQUID loop, which generates the 
magnetic field B(I) at the position of the magnetic particle. 

Equation ([3]) reproduces the results of Refs. [2TII5]. derived for a circular filamentary 
SQUID loop. Moreover, Eq. ([3]) provides a solution of the problem, valid for any 
arbitrary geometry of the superconducting loop, if one can find the normalized magnetic 
field distribution b{f) outside the SQUID loop. For a given b{f) (determined by the 
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Figure 5. (Color online) Calculated spin sensitivity 5 M for SQUID 2 with S$ = 
4/i$ /v / Hz for the detection of the magnetic moment of a small spin particle aligned 
along the x axis in the (x, y) plane of the SQUID loop [c.f. inset in main graph]. Main 

1/2 

graph shows a contour plot of as a function of the position of the particle in the 
(x, z) plane at y = 0. The location of the YBCO bridges is indicated by the black 
rectangles; the rectangles (yellow lines) on top of those indicate the position of the 

1 /2 

Au layer. The lincscan above the main graph shows S^J (x) for z = 85 nm, and the 

1/2 

lincscan right to the main graph shows (z) for x = 0.63 fxm. The location of these 
line scans is indicated by the dashed (red) lines in the main graph. 



SQUID geometry only) and given flux noise, one can use Eq. (jHJ) to easily calculate 
the spin sensitivity Sj/ 2 = S$ 2 /(f)^ for any orientation e M and location of the magnetic 
particle. 

For the geometry of SQUID 2, we calculated the spatial distribution of the 
current density in the SQUID loop and the corresponding 3-dimensional magnetic field 
distribution b(f) outside the SQUID loop with 3D-MLSI |25j. Figure [5] shows the 

1/2 / 

resulting spin sensitivity of SQUID 2 (with S§ = 4/i$ /vHz) for the detection of 
a magnetic particle located in the (x, z) plane (at y — 0) with its magnetic moment 
pointing along the x direction, i.e. e M = e x . I.e., the magnetic moment of the particle is 
aligned parallel to the thin film plane of the SQUID, and perpendicular to the current 
through the GBJs. 

The contour plot of the spin sensitivity shows clear minima right above the 
superconducting bridges straddling the grain boundary. The upper inset shows a 

1 /2 

linescan (x) of the spin sensitivity at a height z = dAu + dy/2 = 85 nm above 
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the ring (i.e. for our SQUID the minimum vertical distance due to the Au layer on 
top of the YBCO film). The lowest value of the spin sensitivity along this linescan 
is 390/ib/vTIz, which could be further improved by reducing the thickness of the Au 
layer. This can be done even without affecting the GBJ properties if the Au layer is 
not removed right above the GBJ. Removing the gold layer (and placing the magnetic 

1 /2 

particle at z — — 25 nm) would improve S M by more than a factor of two down 
to 180 /i^/vTlz, as can be seen in the right inset, which shows the vertical dependence 
S l J 2 {z) at x = 0.63 /xm, i.e. right above the center of the YBCO bridge. Moreover, 
further improvements in S l J 2 are feasible by improving S^ 2 , which is by no means 
optimized for the SQUIDs presented here. For example, our FIB technology allows 
for a reduction in the size of the SQUID loop down to ~ 100 nm and a concomitant 
reduction in SQUID inductance L down to ~ 1 pH. This, in turn, can lead to a significant 

1/2 1/2 

improvement in S$ by at least an order of magnitude, which would bring S/J down 
to « 20/i B /v / Hz. 

5. Conclusions 

In conclusion, we have fabricated YBCO grain boundary junctions and dc SQUIDs by 
FIB patterning with junction widths ranging from 7.8 /im down to 80 nm. Using an 
Au thin film shunt on top of the junctions, we achieved non-hysteretic current-voltage- 
characteristics for operation of YBCO dc SQUIDs at 4.2 K and below. We demonstrated 
that FIB pattering enables the fabrication of deep sub-/im GBJs without degradation 
of critical current densities, and comparable to GBJs with widths above 1 /im. We 
do find a systematic dependence of the resistance times area p of our GBJs, which 
scales approximately with the junction width w as p oc yjw. The origin of this scaling 
could not be resolved and requires further studies. Still, we obtain values of I C R for 
our GBJs around 100 pN for junctions on the 100 nm scale, which is promising for the 
fabrication of sensitive nanoSQUIDs. We demonstrated low-noise performance for such 
devices in the p&o/Vliz range, which still can be improved significantly and which 
makes them promising candidates for applications in magnetic nanoparticle detection 
and measurements at high magnetic fields. The presented solution for calculating 
the spin sensitivity for arbitrary SQUID geometries - as a function of position and 
orientation of the magnetization of small spin particles - provides an important tool 
for the systematic optimization of the spin sensitivity using nanoSQUIDs as sensitive 
devices for direct detection of magnetization switching of small spin particles. 
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